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We present a review on the explanation of the RHIC HBT puzzle by a granular pion- 
emitting source of quark-gluon plasma droplets. The evolution of the droplet is described 
by relativistic hydrodynamics with an equation of state suggested by lattice gauge re- 
sults. The granular source evolution is obtained by superposing all of the evolutions of 
individual droplets. Pions are assumed to be emitted thermally from the droplets at the 
freeze-out configuration characterized by a freeze-out temperature Tf. We find that the 
average particle emission time scales with the initial radius of the droplet. Pions will be 
emitted earlier if the droplet radius is smaller. An earlier emission time will lead to a 
smaller extracted HBT radius Rout, while the extracted HBT radius Rside is determined 
by the scale of the distribution of the droplet centers. However, a collective expansion 
of the droplets can further decrease i?out. As a result, the value of -Rout/Rsidc can be 
close to, or even less than 1 for the granular source of QGP droplets. 



1. Introduction 

HBT (Hanbury-Brown-Twiss) interferometry is an indispensable tool to study the 
space-time structure of the particle-emitting source produced in high energy heavy 
ion collisions Th e exper imental pion HBT measurements at RHIC give the 
ratio of i?out / -Rsidc ~ iHEElTl^ which is much smaller than many earlier theoretical 
expectations. Such a discrepancy between theory and experiment is referred to as 
the RHIC HBT puzzle . On the other hand, hydrodynamical calculations 

give reasonably good descriptions of the elliptic flow, w hich has been c onsidered as 
an evidence for a strongly-coupled quark-gluon plasmgP^ ^ I ^'^ l ^'^ l ^^ l ^^ i The resolu- 
tion of the HBT puzzle is important in finding out why traditional hydrodynamics 
succeed in explaining the elliptic flow but fails in explaining the HBT radii. 

Traditional studies of the hydrodynamics of the evolving fluid assume a sin- 
gle contiguous blob of matter under expansion, with a relatively smooth initial 
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and final density distributions. Initial transverse density fluctuations and hydrody- 
namical instabilities have been neglected but their inclusion may lead to "multi- 
fragmentation" in the form of large scale final-state density fiuctuations and the 
formation of granular droplets. It is useful to explore the consequences of the oc- 
currence of granular droplets. 

Previously we propose a granular model to explain the HBT 

puzzliEdl. We 

would like to review here the important ingredients which enters into the resolu- 
tion of the puzzle. Further suggestions of using single-event HBT interfcrometry to 
search for signatures of the granular source can be found in Refs. '^^ '^^\. 

2. Evolution of a QGP droplet 

Based on the recent results of high-energy heavy-ion collisions at RHIC, the early 
matter produced in the coUisions may be a strongly-coupled QGP (sQGP), which 
has a very high energy density and reaches local thermalization within about 1 
fm/( jl^ | l^ | 14 U^U^. The expansion of the matter after that time may be unsta- 
ble. Many effects, such as the large fluctuations of the initial transverse energy 
density^^ '^'^ the sausage instability^^, and possible phase transitiorPSl^ may lead 
to the fragmentation of the system an d the for mation of many spherical droplets 
due to the surface tension of the QGP I 18 | 19|21 | 

To describe the evolution of a droplet, we use relativistic hydrodynamics where 
the energy momentur n te nsor of a thermalized fluid element in the center-of-mass 
frame of the droplet ij23 

T'^'ix') = [e{x') + pix')]u''{x')u''{x') - p{x')g'"' , (1) 

x' is the space-time coordinate of the fluid element in the center-of-mass frame, e, 
p, and w'^ = 7(1,1;) are the energy density, pressure, and 4-velocity of the element, 
and (7^"^ is the metric tensor. With the local conservation of energy and momentum, 
one can obtain the equations for spherical geometry aJ^ 

dtE + dr[iE+p)v] = -F, (2) 

dtM + dr{Mv+p) = -G, (3) 

where E = T^°, M = T"'\ F = 2v{E+p)/r, G = 2vM/r. 

In the equations of motion ([2]) and ([3]) there are three unknown functions e, p, v. 
In order to obtain the solution of the equations of motion, we need an equation of 
state which gives a relation p{e) between p and e ■^^J . At RHIC energy, the system 
undergoes a transition from the QGP phase to hadronic phase. As the net baryon 
density in the central rapidity region is much smaller than the energy density of 
the produced matter (here presumed to be QGP), the baryon density of the system 
in the center rapidity region can be neglected. Lattice gaug e results s uggest the 
entropy density of the system as a function of temperature aJSIMHHSni 



siT)_[T\'\^^d^^fT-T. 



dg + dH \ AT 



(4) 
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where Sc is the entropy density at the transition temperature Tc, dq and dn are 
the degrees of freedom in the QGP phase and the hadronic phase, and AT is the 
width of the transition. The thermodynamical relations among p, e, and s in this 
case are 

—sdT + dp = 0, e = Ts—p. (5) 

From these thermodynamical relations and Eq. ([4]), we can obtain the equation of 
state p{e). 

Using the HLLE schemePll^ and Sod's 
operator splitting method"^, on e can ob tain 
the solution of Eqs. Q and ([SlP^SIIZ!^ af- 
ter knowing the equation of state and initial 
conditions. We assume that the droplet has 
a uniform initial energy density eg within a 
sphere with radius r^, and has a zero ini- 
tial velocity in its center-of-mass frame. Figs. 
1(a) and (b) show the temperature profiles 
and isotherms for the droplet. In our calcu- 
lations, we take the parameters of the equa- 
tion of state as dq = 37, dn = 3, Tc = 165 
MeV, and AT = 0.05 T^, and take the initial 
energy density eg = 3.75TcSc, which is about 
two times of the density of quark matter at 

3. Granular Source of QGP droplets 

If we assume that the final pions are emitted 
from the droplet at the freeze-out configura- 
tion characterized by a freeze-out tempera- 
ture Ty, we can see from figure 1(b) that the the average particle emission time 
scales with the initial radius of the droplet r^. In HBT interferometry, the radius 
i?sidc is related to the spatial size of the particle-emitting source and the radius 
-^out is related not only to the source spatial size but also to the lifetime of the 
gQ^j.(,(J213134135| ^ jQj^g lifetime of the source will lead to a large Ron^^^^^- From 
the hydrodynamical solution in figure 1(b), both the average freeze-out time and 
freeze-out radial distance increase with Td for a single droplet source. As a conse- 
quence, -Rout /-Rside is inscnsitivJI^ to the values r^. The value of i?out /-Rside forthe 
single droplet sourcJI^ is about 3 Pj , much larger than the observed valueJ3EISEl_ 
The RHIC HBT puzzle of i?out / -Rsidc ^ 1 suggests that the pion emitting time 
may be very short. In order to explain the HBT puzzle, we consider a granular source 
of Nd QGP droplets distributed in a distribution D{Xd) as illustrated in Fig. 2(a). 
We assume that the evolution of the granular source is simply the superposition of 



1.6 
1.2 
0.8 
0.4 

6 



4 



— 1 — I — I — I — I — I — I — 
U (a) 



tn \ 



I I I I I ' 



o.sr^ 



• ^o.65r,. 



(b) 



0.95 r, 



•, 1 

: ( 

•/ 
•/ 
•/ 

■y 

_] I I I L_ 



r/r^ 

Fig. 1. (a) Temperature profile and 
(b) isotherms for the droplet. Here, 
t„ — inXrd and A — 0.99. 
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the evolutions of individual droplets with the same initial conditions. As the average 
freeze-out time is proportional to the initial radius of the droplet, the freeze-out time 
and i?out decreases if the initial radius of the droplet decreases. On the other hand, 
i?sido increases if the width of the droplet spatial distribution D{Xd) increases. 
A variation of the droplet size and the width of droplet spatial distribution can 
result in i?out nearly equal to i?sidG l^'*']- Furthermore, if the granular source has a 
collective expa nsion the HBT will measure the size of the region of the ellipse in 
figure 2(b)^^. In this case the value of i?out/-Rsidc decreases even further and may 
be smaller than unityl^. 




(a) (b) 

Fig. 2. (a) Static Granular source of the QGP droplets, (b) Expanding granular source of the QGP 
droplets. 



We assume that the droplets are initially distributed in a short cylinder of 
length 2TZz along the beam direction (z direction) with an initial transverse spatial 
distribution up to a radius TZt, 

[1 - exp(-pVA7^?)]0(7^, - p) 9{TZ. - \z\) , (6) 



271 p dp dz 



where p = x'^ + and z are the coordinates of the center of a droplet, and 
ATZt describes a shell-type radial distribution. Because of the early thermalization 
and the anisotropic pressure gradient, the droplets will acquire anisotropic initial 
velocities. We therefore consider a granular source of the the QGP droplets with an 
anisotropic velocity distribution and the velocity of a droplet depends on the initial 
coordinates of the droplet center, (?'i,r2,r3) = (x,y,z), in the form 



f3^ = a, sign(r,) j . « = 1, 2, and 3, (7) 

where Oi describes the magnitude of the anisotropic expansion, sign(ri) denotes 
the sign of r^, and bi (6^, bz) are the exponential power parameters that describe 
the variation of the velocities with r^. In our calculations the velocity parameters 
are taken as = 0.415, ay — 0.315, Qz = 0.850, b^ = by = 0.42, and bz = 0.03, 
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by comparing with the experimental data of pion transverse momentum spectrai^^l 
and elliptic flo-wEIlin ^/snn = 200 GeV Au + Au collisions at RHIC. We take the 
freeze-out temperature as Tf = 0.95Tc in our calculations. 

4. HBT Results of the Granular Source of QGP Droplets 

The two-particle Bose- Einstein correlation function is defined as the ratio of the two- 
particle momentum distribution P(pi,p2) relative to the the product of the single- 
particle momentum distribution P(pi)P{p2). For a chaotic pion-emitting source, 
P{pi) {i = 1,2), and P{pi,p2) can be expressed aJU 

P(p,)=5]a2(p„X,), P(pi,P2)= J2 \Hpi,P2;X,,X2)\\ (8) 

Xi Xi,X2 

where A{pi,Xi) is the magnitude of the amplitude for emitting a pion with 4- 
momentum pi = (jj^, Ei) in the laboratory frame at Xi and is given by the Bose- 
Einstein distribution with freeze-out temperature Tf in the local rest frame of 
the source point. ^{pi,p2; Xi, X2) is the two-pion wave function. Neglecting the 
absorption of the emitted pious by other droplets, ^{pi,p2; Xi, X2) is simply 

$(pi,P2;^i,^2) = ^[A{pi,Xi)A{p2,X2)e''"-'''+''''-''-' + (X, ^ X2)] . (9) 

Using the components of "out" , "side" , and "long"l2ll25| Qf i^j^g relative momentum 
of the two pions, q = |pi — P2I, as variables, we can construct the correlation 
function C(gout, gsido, 9iong) from P{pi,p2) and P{pi)P{p2) by picking pion pairs 
from the granular source and summing over pi and P2 for each (gout, Qside, 9iong) 
The HBT radii i?out, ^sido, and i?iong can then be extracted by fitting the 
calculated correlation function C(gout, ?side) Qiong) with the following parametrized 
correlation function 

C(gout, gsidc, giong) = 1 + A e-'?°-«=--'?a-«'--«?-.<". . (10) 

Figure 3 shows the theoretical two-pion correlation functions for the granular 
source. The top figures give the results for the average pion transverse momentum 
of a pion pair, Kt^ less than 400 MeV/c, and the bottom figures for Kt > 400 
McV/c. In our calculations the size parameters of the granular sources are taken 
to be TZt ~ 8.8 fm, AT?.* — 3.5 fm, TZz — 7.0 fm, and ra — 1.3 fm. The number 
of droplet is taken to be 40 in our calculations. The left panels in Fig. 4 give 
the extracted two-pion HBT radii for the granular source as a function of Kt- 
The symbols of circle and down-triangle are for /SJZt — 0.35 fm and A7?,t = fm, 
respectively. The experimental PHENIX resultsP, and STAR resultJ^ are shown 
on the right panels. The curve gives the theoretical results for MZt = 0.35 fm. 
For our theoretical HBT calculations, we use a cut for particle pseudo-rapidity 
region \ri\ < 0.35, the same as in the PHENIX experimentJ^. We find that if we 
increase the parameter r^, the HBT radii -Rout and -Riong will increase. And if we 
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Fig. 3. Two-pion correlation functions for granular source. 

increase the parameter A7?.t, the variation of HBT radii i?out and i?sidc with Kt 
win become steep. As can be inferred from Fig. 4, the HBT results of granular 
source for A7?,( = 3.5 fm agree quite weh with experimental data, although the 
case with MZt = also give almost as good an agreement. 



5. Conclusions and Discussion 

The expansion of the dense matter (sQGP) produced in high-energy heavy-ion 
collisions at RHIC may be unstable. The initial transverse density may also be 
highly fluctuation. The unstable expansion and fluctuating initial transverse density 
may lead to a fragmentation of the system and the formation of a granular source of 
QGP droplets. Although a granular structure was suggested earlier as the signature 
of first-order phase transitionJ^, the occurrence of granular structure may not be 
limited to the occurrence of first-order phase transitions. There ar e addi tional effects 
which may lead to the dynamical formation of granular 

dropletJlSEII 

For a granular source of the droplets, the average particle emission time scales 
with the initial radius of the droplet. Pions will be emitted earlier if the radius is 
smaller. An earlier emission time will lead to a smaller extracted HBT radius i?out, 
while the extracted HBT radius i?side is determined by the scale of the distribution 
of the droplets. A collective expansion of the droplets can further decrease i?out 
and the value of i?out/-Rsidc can be close to, or even less than 1 for an expanding 
granular source of QGP droplets. 

In an event-mixing experimental analysis, some signals of the granular source 
such as the correlation function fluctuations may likely be suppressed after averag- 
ing over many eventJ^^. However, the property of an earlier average emission time 
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Fig. 4. Two-pion HBT radii obtained by the PHENIX Collaboration^ and the STAR 
Collaboration'^ compared with the theoretical results calculated in the granular droplet model 
with ATZt = and 3.5 fm. 



of the granular source will remain as it is common to all events and not averaged 
out by event mixing. This property leads to a smaller i?out/^sidc in mixed-event 
HBT analysis. 

In conclusion, a granular source model can explain the RHIC HBT puzzle and 
reproduce the data of pion HBT rad ii, as well as the data of pion transverse mo- 
mentum spectra and elliptic flo'vJl^ in ^/snn = 200 GeV Au -I- Au collisions at 
RHIC. It is of great interest to find direct evidence s of the granula r structure and 
to study the mechanics of granular source 

formatioi|2IIlI19l2Ql2T]_ 
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